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Abstract 

We study the statistical properties of complex networks constructed from time series of energy dissipation rates in 
three-dimensional fully developed turbulence using the visibility algorithm. The degree distribution is found to have 
a power-law tail with the tail exponent a — 3.0. The exponential relation between the number of the boxes Nb and the 
box size lg based on the edge-covering box-counting method illustrates that the network is not self-similar, which is 
also confirmed by the hub-hub attraction according to the visibility algorithm. In addition, it is found that the skeleton 
of the visibility network exhibits excellent allometric scaling with the scaling exponent rj = 1.163 + 0.005. 
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scaling 
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1. Introduction 

Turbulence is one of the most challenging fields in physics, and diverse methods have been adopted to investigate 
the statistical properties of turbulent flows QjJ, . For instance, the multifractal nature of turbulent flows has been ex- 
tensively studied based on the structure function analysis of velocity time series ^ and the partition function analysis 
of energy dissipation rates J3,|5t]. Due to the extreme complexity of turbulence, phenomenological investigation of 
experimental data plays a crucial role in order to gain a better understanding of turbulence, which is fundamental for 
the construction of models and theories. In this work, we provide a first attempt to study turbulent signals from the 
complex network perspective, hoping that the ideas, tools, and theories from the flourishing field of complex networks 
can stimulate the study of turbulence from an alternative point of view. It has been shown that the idea of complex 
network analysis is able to identify flow patterns and nonlinear dynamics of gas-liquid two-phase flows Jf^O]. 

Recently, several methods that convert time series into networks have emerged. Zhang et al. introduced a method 
to deal with the pseudoperiodic time series and found that the structure of the corresponding network depended on the 
dynamics of the series fllgt]. In this method, the nodes correspond directly to cycles in the time series, and the edges 
are determined by the strength of temporal correlation between cycles. This method for pseudoperiodic time series 
can also be generalized to other time series, where a node is defined by a sub-series of fixed length as an alternative of 
a cycle, which has been applied to stock prices ifioll . Xu et al. proposed another method, which embeds the time series 
into an appropriate phase space, takes each phase space point as a node in the network, and connects each point with 
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its four nearest neighbors to form a complex network ll 111 . There are also other mapping methods from time series to 
networks based on n-tuples of fluctuations 1 12 - 14 ] or recurrence plots 1 15 , . 

Lacasa et al. proposed the visibility graph algorithm, which can map all types of time series into networks jV, 
In the visibility algorithm, the nodes correspond to the data points of the time series in the same order, and an edge 
is assigned to connect two nodes if they can "see" each other. It is found that the degree distributions of visibility 
graphs converted from self-similar time series have power-law tails |fl7ll and the tail exponents depend linearly on the 
Hurst index of the original time series 1 1_8 , lfj, in which the multifractal nature of the time series has a negligible 
influence I 1 8 |j . The visibility algorithm can be simplified to a variant termed the horizontal visibility algorithm, which 
is solvable H20fl . The visibility algorithm has been diversely used to investigate stock market indices Il8[|2lh. human 
strike intervals Jl9ll . the occurrence of hurricanes in the United States 1221. and foreign exchange rates ll23h . 

In this paper, we study the statistical properties of complex networks constructed from the energy dissipation rate 
time series using the visibility algorithm. In Section|2] we describe the turbulence data set. Section[3]briefly describes 
the visibility graph algorithm. The main results on the degree distribution, the nonfractality and the allometric scaling 
of the constructed visibility are given in Section^] We summarize in Section[5] 



2. Description of the data set 

The velocity data set from three-dimensional fully developed turbulence has been collected at the SI ONERA 
wind tunnel by the Grenoble group from LEGI 0]. The mean velocity of the flow is approximately (v) = 20m/s 
(compressive effects are thus negligible). The root-mean-square velocity fluctuations is v rms = 1.7m/s, leading to a 
turbulence intensity equal to / = v rms /(v) = 0.0826. This is sufficiently small to use Taylor's frozen flow hypothesis. 
The integral scale is approximately 4m but is difficult to estimate precisely as the turbulent flow is neither isotropic 
nor homogeneous at these large scales. The Kolmogorov microscale 77 is given by J^t] 



n = 



2/ \2 n 1 / 4 
v <v}~ 

= 0.195 mm, (1) 



15((dv/dt) 2 ) 

where v = 1.5 x 10 _5 m 2 s _1 is the kinematic viscosity of air, and dv/dt is evaluated by its discrete approximation with 
a time step increment dt = 3.5466 x 10~ 5 s corresponding to the spatial resolution 5 r = 0.72mm divided by (v), which 
is used to transform the data from time to space by applying Taylor's frozen flow hypothesis. The Taylor scale is given 
by J! 

, (v)v rms 

The Taylor scale is thus about 85 times the Kolmogorov scale. The Taylor-scale Reynolds number is 

Re,! = Av^/v = 2000. (3) 

This number is actually not constant along the whole data set and fluctuates by about 20%. We have checked that 
the standard scaling laws reported in the literature are recovered with this time series. In particular, we have verified 
the validity of the power-law scaling E(k) ~ with an exponent /? very close to 5/3 over a range more than two 
orders of magnitude, similar to Fig. 5.4 of [Q]] provided by Gagne and Marchand on a similar data set from the same 
experimental group. 

Then, we obtained the kinetic energy dissipation rate data using Taylor's frozen flow hypothesis which replaces 
a spatial variation of the fluid velocity by a temporal variation measured at a fixed location, and the rate of kinetic 
energy dissipation at position i is 

ei~[{v M -Vi)l8 t ? (4) 

where 6c is the resolution (translated in spatial scale) of the measurements. Fig.Q]shows a randomly selected segment 
of 4000 values of the energy dissipation rate time series. The energy dissipation rate time series exhibits multifractal 
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Figure 1 : A segment of 4000 data points of the energy dissipation rate time series. 



3. Construction of visibility graph 

Here, we give a brief introduction of the visibility graph algorithm. For a given time series (?, e), where t is the 
time variable and e is the value of energy dissipation rate, a visibility line exists between two points (t a , e a ) and (tb, €b), 
if any other data (t c , e c ) placed between them fulfills: 

e c < e h + (e a - <%)^ — - . (5) 

tb ~ ta 

Fig. |2] illustrates the scheme of the visibility algorithm, where eight points of the energy dissipation rate series 
are plotted as an example. The height of the vertical lines are the values of energy dissipation rates. The cycles and 
the dashed lines constitute the visibility graph. The nodes correspond to series data in the same order and an edge 
connects two nodes if one can be seen from the top of the other. In other words, two nodes are connected if there is 
visibility between them. Take points 1 and 4 of Fig. [2]as the example to explain the concept of "visibility". Between 
1 and 4, there are two points 2 and 3 which are all under the dashed line from point 1 to point 4, and visibility exists 
between nodes 1 and 4. 
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Figure 2: The visibility graph of a short time series. 

The networks mapped according to the visibility graph algorithm are connected and undirected. Different time 
series convert into different networks. It is easy to check that for the time series with monotonously decreasing slopes 
of all the adjacent points, a chain-like graph will be obtained. According to this algorithm, the dynamics of the time 
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series is conserved in the graph topology. For instance, Lacasa et al. found that periodic time series convert into 
regular graphs, random time series into random graphs, and fractal series into scale-free graphs lfl7ll . 



4. Statistical properties of the visibility graph 

4.1. Degree distribution 

Fig.[3]illustrates the degree distribution p(k) of the visibility graph mapped from the energy dissipation rate series 
in three-dimensional fully developed turbulence with 10 6 data sets. It is evident that the degree distribution has a 
power-law tail 

P(k) ~ k~ a , (6) 

where the exponent is estimated to be a — 3.0. This power-law distribution is obtained empirically using the method 
proposed by Clauset, Shalizi and Newman |2(J. The power-law exponent is determined by maximum likelihood 
estimation based on the Kolmogorov-Smirnov statistic 12611 and the power-law distribution is statistically significant. 
It means that the visibility graph under investigation is scale-free. The power-law degree distribution has been found 
in the visibility graphs mapped from fractal series such as the Brownian motion and Conway series lfl7ll . Fig.[3]shows 
that multifractal data sets of the energy dissipation rate in turbulence also converts into scale-free networks and the 
exponent of the degree distribution of the network may depend on the correlation of the series. 




Figure 3: The degree distribution of the visibility graph converted from the time series of the energy dissipation rate series. The line is the best 
power-law fit of the data. 



The power-law degree distribution of the visibility graph has implications for the fractality of the original time 
series Il7l 4l9]l . This is consistent with the well known result that the energy dissipation rate series is statistically 
self-similar [1]. For nonstationary time series, such as stock price indexes, fractional Brownian motions, Conway 
series, and strike intervals, there is a linear dependence of the Hurst index on the tail exponent 1 18, lfj. However, no 



quantitative relation between the Hurst index of stationary time series and the tail exponent of the degree distribution 
has been established. Indeed, if we apply rudely the relation in Refs. JI8, 19], we will come to the conclusion that the 
energy dissipation rate series is anticorrelated since the resulting Hurst index is less than 1/2, which is consistent with 
the empirical facts 1271- 



4.2. The visibility graph is not self-similar 

The relation between scale-free networks and self-similar networks has been discussed recently 1 28 -30j], and 
fractal networks seem to exhibit scale-free features while scale-free networks are not always self-similar. The fractal 
nature of the network can be revealed by the well-known box-counting method. Either the node-covering [28] or the 
edge-covering [31] box-counting method needs to calculate the minimum number Nb of boxes with size Ib that can 
cover all the nodes or edges of the network, where Is - 1 is the maximum distance of all possible pairs of nodes in each 
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box in the node-covering box-counting method while Ib represents the largest distance in the edge-covering method. 
If the network is self-similar, the minimum number Nb of covering boxes scales with respect to box size Ib as 



(7) 



where D is the fractal dimension of the network. In order to obtain the minimum number of the covering box, the 
simulated annealing algorithm is implemented i3lll and the computation may be too enormous for the network mapped 
from the whole series. In this work, the graph mapped from the 1000 continuous points of the turbulence series are 
considered to study the fractal feature. Fig. [4] shows the results of the edge-covering box-counting analysis of the 
visibility graph. As illustrated in Fig. [4] the number of boxes Nb decreases exponentially with the size of the boxes 
Ib, rather than a power law. Thus, the visibility graph in our work is not fractal. 




Figure 4: Edge covering box-counting analysis of the visibility graph. Ng and Ib are in an exponential relation, and the network is non-fractal. 



According to Song et al. 1281 12911 . the network with hubs (the high-degree nodes) attraction will not be self- 
similar. Whether the interaction is attraction or repulsion among the hubs of the network can be determined from 
the correlation between the degrees of different nodes. The degree correlation can be quantified by a correlation 
coefficient r 11321 13311 : 

ihk 2 ) - {h){k 2 ) 

r = -, (8) 



where cr^ is the variance of all the node degrees, and k\, k 2 are the degrees of the nodes at the ends of all the edges. 

It is clear that the degree correlations of adjacent nodes have strong effects on the structure of the network. Net- 
works with positive correlation (r > 0) tend to have a core -periphery structure, and the hubs are attracted to each other 
with large probability, while the hubs of a network with negative correlation (r < 0) are repulsive. In our case, we 
obtain that r = 0.1204 for the network mapped from the whole time series. The positive correlation of the whole net- 
work indicates that the graph mapped from the energy dissipation rate series is hub-attractive and thus not self-similar. 
According to the visibility graph method, it is easy to find that large values of the energy dissipation rates correspond 
to the hubs in the network with large probability, and the points with large values can "see" each other easily which 
means that edges exist between the hubs. 

More precisely, the correlation coefficient r in Eq. © is not a good indicator to determine the fractality of the 
network. As is explained in Ref. l29tl . for the property of a higher degree of hub repulsion to be the hallmark of 
fractality, "the anticorrelation has to appear not only in the original network, but also in the renormalized networks 
at all length scales." What is important is not the degree correlation but the scaling of the degree correlation l34ll . In 
addition, some measures of anticorrelation, such as the Pearson coefficient r, cannot capture the difference between a 
fractal and a non-fractal network, since r is not invariant under renormalization J29 , 34 ] . However, positive correlation 
of degrees using the Pearson coefficient is sufficient to confirm that the network is not fractal. 
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4.3. Allometric scaling of the visibility graph 

The allometric scaling was first proposed in biology 1 35, 36|, and it turns out that the various scalings are charac- 
terized by the allometric equation: 

Y ~ M b (9) 

where Y is any physiological, morphological or ecological variable that appears to be correlated with size or body 
mass M, and b is the scaling exponent. Allometric scaling laws have been uncovered in many networks, such as 
organ metabolic networks 1 37], food webs lf38h . world trade/investment webs 1 39[ 4(J, and so on |41]. The allometric 



scaling law is proposed in a network based on spanning trees. With one root and plenty of branches and leaves, a 
tree can be considered as directed from root to leaves. For each node i in the tree, we calculate A,- of node ; that the 
summation of value A of its branches along with node i itself 



I> +1 



(10) 



and the quantity C, of node ; is the sum of the C values of its branches and A, 

C i = Y A Cj+A u (11) 

,/' 

where j is the son node of ;. 

Imagining that there is a flow from the root to the leaves of the tree, the value A stands for the size of node while 
the value C represents the transport rate through the node. According to Eq. (01, if the property of allometric scaling 
exists in the tree, the rate C, and the size A, should depend on a power-law relation: 

C~A n . (12) 

The exponent 77 is a very important parameter of the tree structure. The transportation is more efficient for smaller 
77 values. The chain-like trees and star-like trees are the two extremes of the spanning trees, and the value of the 
allometric exponent 77 is also between that of the two trees. It is easy to obtain that 77 = 2 for chain-like trees and 
77 - 1 + for star-like trees. It follows that 1 < 77 < 2. We note that not all trees exhibit allometric scaling properties 



H39U42U . and different spanning trees for the same network may exhibit different allometric features. 

For a network, there are various spanning trees with very different structures. In this work, we consider the 
skeleton 1 3(1 43 ] of the network. The skeleton i s a p articular spanning tree which is formed by the edges with largest 



betweenness centrality of the original network 143ll . The skeleton of the scale-free network is also scale-free with a 
different exponent a, and if the original network is self-similar, the skeleton also exhibits self-similarity. To extract 
the skeleton of a network, we first determine the shortest paths between all pairs of nodes, and the number of times an 
edge appears in the shortest paths is defined as the weight of the edge, 



(E) = Y j c(i,j,E), (13) 



>J 

where c(i, j, E) — 1 if the shortest path from node ; to j through edge E and c(i, j, E) = otherwise. The node with the 
largest degree is considered as the root of the spanning tree. As a first step, we select an edge with the largest weight 
and add it to the tree if its inclusion will not form any loops. Once all the nodes of the original network are added into 
the tree, the skeleton is obtained. 

We randomly select 100 segments of the time series, each segment having 1000 data points. The 100 sub-series 
are converted into visibility graphs and their skeleton spanning trees are determined accordingly. Fig. [5] shows the 
comparison of the scaling behavior between the skeleton and a random spanning tree. The plot of skeleton is shifted 
upward for clarity and the lines are the best fits of the data. The perfect power-law relations are observed between C 
and A of all the considered segments, and the line is the best fit of In C against In A. The leaves of the trees (A = 1, 
C — 1 ) are excluded from the fitting i38ll . The exponents of all the skeletons of different segments are around 1 . 1 6 with 
slight fluctuations, indicating that the allometric scaling property of the visibility graph skeleton remains unchanged at 
different parts of the turbulence time series. We find that 77 = 1.163+0.005. We also observe that the data points of the 
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random spanning tree are more dispersive in the scatter plot than those of the skeleton. It is interesting to note that this 
allometric scaling exponent is quite universal in many different networks 11381 14211 . However, the physical implication 
and the relation to the dynamics of the turbulence signal of this universal topological structure of the visibility graph 
are unclear. 
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Figure 5: Comparing the allometric scaling behavior of the skeleton and a random spanning tree of the visibility graph mapped from the same 
series. The lines are the best power-law fits of the data. The data for the skeleton have been shifted upward by two orders of magnitude for clarity. 



5. Summary 

The visibility graph algorithm is a new method for time series analysis, which converts time series into complex 
networks. In this way, the time series can be studied from a complex network perspective. The visibility graphs inherit 
many properties of the associated time series. We have investigated the visibility graphs constructed from the time 
series of energy dissipation rate in turbulence. We found that the visibility graph has a power-law tail in its degree 
distribution and is thus scale free, which corresponds to the self-similarity of the energy dissipation rate series. The 
edge-covering box-counting method shows that the number of boxes decays as an exponential function with respect 
to the box size, which indicates that the visibility graph is not fractal. The positive degree correlation coefficient r and 
the visibility algorithm for the turbulence data show hub-hub attraction (rather than hub-hub repulsion) in the network, 
which explains why there is no self-similarity. In addition, the allometric scaling properties of the visibility graph are 
studied. We selected skeleton as the particular spanning tree of the network, and the results show that the skeleton 
exhibits excellent allometric scaling. The networks mapped from various segments of the series of the turbulence have 
the same scaling exponent 77 = 1.163 + 0.005. However, the physical implication of the allometric scaling is not clear. 
We hope that our analysis will stimulate the study of turbulent signals from the perspective of complex networks. 
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